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Abstract 

A new relativistic Hartree-Fock approach with density-dependent a, to, p and 7r meson-nucleon couplings for 
finite nuclei and nuclear matter is presented. Good description for finite nuclei and nuclear matter is achieved 
with a number of adjustable parameters comparable to that of the relativistic mean field approach. With the 
Fock terms, the contribution of the 7r-meson is included and the description for the nucleon effective mass and 
its isospin and energy dependence is improved. 

PACS numbers: 21.60.-n, 21.30.Fe, 21.60.Jz, 21.10.Dr, 24.10.Cn, 24.10.Jv 


The relativistic mean field (RMF) theory 
description of numerous nuclear phenomena 



31 has received much attention due to its successful 
. In its most widely employed versions, 


a? imi? LLii 


i.e., with either self-coupling interactions or density-dependent meson-nucleon couplings, the RMF 
theory with a limited number of parameters can describe very well a very large amount of data: 
saturation pro per ties of nuclear matter |icj], nuclear binding energies and radii, the isotopic shifts in 


the Pb-region 


iperti 

fi. 


It gives a natural description of the nuclear spin-orbit potential 




ial Q, 


and explains 


the origin of the pseudospin symmetry [131.1141] and spin symmetry of the anti-nucleon spectrum [15| as 


a relativistic symmetry |15I. lid. 1171 lla|. In spite of these success, there are still a number of questions 
needed to be answered in the RMF theory: the contributions due to the exchange (Fock) terms and 
the pseudo-scalar 7r-meson. 

There exist attempts to include the exchange terms in the relativistic description of nuclear matter 
and finite nuclei. The earlier relativistic Hartree-Fock (RHF) method led to underbound nuclei due to 


the missing of the meson self-interactions 12] • Further developments were made by taking into account 


approximately the nonlinear self-couplings of the c-field m or by introducing the products of six and 
eight nucleon spinors in the zero-range limit 0I| . Although some improvements were obtained, the 
RHF method is still not comparable with the RMF theory in the quantitative description of nuclear 
systems. The relativistic point coupling model has been used to investigate nuclei systems Q] and the 
consequences of Fierz transformations acting upon the contact interactions for nucleon fields occurring 
in relativistic point coupling models has been investigated in Hartree approximation, which yield the 
same models but in Hartree-Fock approximation instead QQ. It has been suggested that the 
Hartree-Fock approximation may constitute a physically more realistic framework for power counting 
and QCD scaling than the Hartree approximation. 

In this work, a new RHF approach which contains density-dependent meson-nucleon couplings is 
developed. With a number of adjustable parameters comparable to that of RMF Lagrangians, this 
density-dependent RHF (DDRHF) theory can give a good description of nuclear systems without 
dropping the Fock terms. Furthermore, important features like the behavior of neutron and proton ef¬ 


fective masses 25i| can be interpreted well in DDRHF in comparison with the results of non-relativistic 

n nn 

Brueckner-Hartree-Fock (BHF) [26j] and Dirac-Brueckner-Hartree-Fock (DBHF) calculations 127. i:2a]. 

The most important parts of the nuclear force are the short-range repulsive and medium-range 
attractive components. In analogy with the strong interaction in free space which is described by meson 
exchanges, it is convenient to represent the strong interaction in nuclear medium by the exchange of 
effective isoscalar and isovector mesons. The description of nucleon and meson degrees of freedom 
has to rely ultimately on the relativistic quantum held approach. According to this spirit, we start 
from an effective Lagrangian density 2z? constructed with the degrees of freedom associated with the 
nucleon held(^), two isoscalar meson fields (a and u), two isovector meson fields (it and p) and the 
photon held (A). The parameters of the model are the effective meson masses and meson-nucleon 











couplings. 


With the general Legendre transformation 


= T 00 = - Jgf, 

dfa 


one can obtain the effective Hamiltonian from the Lagrangian density 2z? as 


(1) 


=ip {—i'Y ■ V + M) if; 

+ i / d 4 X 2 £ ip(xi)i/;(x2)TiDi(xi, x 2 )ip(x2)ip{xi), ^ 

i=(T,LO, 

P,7T,A 

where D{(xi,x 2 ) represent the corresponding meson propagators, and the interaction vertices T* are 
defined as, 


T ff (l,2) = - g ff (l)g ff (2), (3a) 

r w (l, 2) = + ^(l) 7M (l)^(2) 7 /i (2), (3b) 

F P (1, 2) = + g p {l)^{l)r{l) • 5p (2) 7 ^(2)f(2), (3c) 

r„(l,2) = - [^757^5^- [^f 7 57^] 2 , (3d) 

La(1, 2) ee + j [ 7m (1 - rs)], [7^(1 - r 3 )] 2 . (3e) 


Following the experience and success in DDRMF |8|, |2H l30l . I31l. I32l . 1331 . |3J|, the meson-nucleon 
couplings (ja■ (Juj ■ g P and /, r are taken as functions of the baryonic density pb- For a- and w-meson, the 
density-dependence of the couplings g a and g w are chosen as 


giiPb) = gi(po)fi(0, 


(4) 


where i = a, u, and 


fi{0 = ai 


1 + &i(£ + d*) 2 
1 + Cj(^ + dj) 2 


(5) 


is a function of £ = pb/po, and po denotes the baryonic saturation density of nuclear matter. In 
addition, five constraint conditions /)(1) = 1,/"(1) = f”{ 1), and /"(0) = 0 are introduced to reduce 
the number of free parameters. For simplicity, the exponential density-dependence is adopted for f n 
as well as g p [35]: 


9 P {Pb) =g P { 0)e apC , 

(6a) 

fn(Pb) =fn(0)e~ a7r ^. 

(6b) 


The coupling constants g p ( 0) and f 7r ( 0) are fixed to their values in free space. One reason to fix the 
coupling constants g p ( 0) and f n ( 0) to their values in free space is just to reduce the number of free 
parameters and another reason is that the inclusion of Fock terms allows such choice. There are in 










total 8 free parameters, i.e., m a , g a (po ), g^iPo), a pi a m and three others from the density-dependence 
of g a and g u . A new parametrization called PK01 is found (see Table EJ by fitting the masses of the 
nuclei 16 0, 40 Ca, 48 Ca, 56 Ni, 68 Ni, 90 Zr, 116 Sn, 132 Sn, 182 Pb, 194 Pb, 208 Pb and 214 Pb, and the values of 
the baryonic saturation density po , the compression modulus K and the symmetry energy J of nuclear 
matter at the saturation point. 

It should be emphasized that here the effective interaction PKOl is obtained by fitting the empirical 
properties of nuclei and the nuclear matter at the saturation point. In Table I one finds a p = 0.076 
and a n = 1.232. This means that g p (l)/g p (0) =0.93 and /7r(l)//7r(0) =0.36, i.e., the contribution from 
the pion is reduced heavily than that in free space. In fact, the effect of pion has been taken into 
account effectively via the other mesons. It will be refined in the future if more information is used 
to constrain the density dependence of the effective interaction in the medium. 


TABLE I: The effective interaction PKOl for DDRHF with M = 938.9MeV, m u = 783.0MeV, m p = 769.0MeV 
, = 138.0MeV. 


m a 

525.7691 

CL<j 

1.3845 

a u 

1.4033 

9a 

8.8332 

ba 

1.5132 

bu 

2.0087 

9uj 

10.7299 

C-cr 

2.2966 

C(jj 

3.0467 

9p( °) 

2.6290 

d(j 

0.3810 

du 

0.3308 

U 0) 

1.0000 

a p 

0.0768 

Cl-jr 

1.2320 


The PKOl parameter set gives the following nuclear matter bulk properties: compression modulus 
K = 250.24MeV, symmetry energy J = 34.37MeV, binding energy per particle E/A = —15.996MeV, 
saturation baryonic density po = 0.1520fm -3 . 

For finite nuclei, the self-consistent Dirac equations are solved in coordinate space with techniques 
similar to those used in RMF 2^,1^]. The non-local exchange (Fock) potentials are treated exactly as 


in Ref. 


Q. 


Calculations are carried out for a set of selected nuclei (S.N.), i.e., 16 0, 40 Ca, 48 Ca, 56 Ni, 


58 Ni, 68 Ni, 90 Zr, 112 Sn, 116 Sn, 124 Sn, 132 Sn, 182 Pb, 194 Pb, 204 Pb, 208 Pb, 214 Pb, and 210 Po, as well as the 
Sn and Pb isotopic chains. For the open shell nuclei, the pairing correlations are treated by the BCS 


method with a density-dependent delta force J38j. A detailed comparison with the predictions of some 


typical RMF parameterizations: PK1 j^j], PKDD j^j], NL3 3£] and DD-ME1 |33| are summarized 
in Table ED where the root mean square (rms) deviations from the data are shown. As one can see in 
Table ED the DDRHF approach with PKOl provides a good quantitative description of finite nuclei, 
sometimes better than the RMF approach. It should be emphasized that this is the first time for the 
RHF approach to provide such a good quantitative description for the finite nuclei and nuclear matter. 

From previous discussion, one can find that good description of nuclear systems comparable to that 
of RMF can be obtained without dropping the Fock terms. Taking 208 Pb as an example, shown in 
Fig. HI are the neutron energy densities from Hartree and Fock terms in different meson channels in 
















TABLE II: The rms deviations A from the data for the RHF calculations with PKOl in comparison with that 
of RMF with PK1, PKDD, NL3 and DD-ME1. The rows from two to ten are respectively: the binding energies 
Eh of the selected nuclei (S.N.) and the even-even nuclei in Pb and Sn chains; two-neutron separation energies 
S‘ 2 n of Pb and Sn isotopes; charge radii r c of the S.N. and Pb isotopes; isotope shifts (I.S.) of Pb isotopes; 
spin-orbit (S.O.) splittings of doubly magic nuclei. 




PKOl 

PK1 

PKDD 

NL3 

DD-ME1 


S.N. 

1.6177 

1.8825 

2.3620 

2.2506 

2.7561 

A E b 

Pb 

1.8995 

2.0336 

2.7007 

2.0021 

2.1491 


Sn 

1.2665 

1.9552 

2.4567 

1.6551 

0.9168 

A S 2 „ 

pb 

0.6831 

0.9192 

1.3139 

0.9359 

1.2191 

Sn 

0.6813 

0.7762 

1.0629 

0.8463 

0.7646 


S.N. 

0.0269 

0.0204 

0.0188 

0.0177 

0.0163 

A r c 

Pb 

0.0056 

0.0061 

0.0060 

0.0143 

0.0150 

Ai.s. 

Pb 

0.0760 

0.0784 

0.0784 

0.0679 

0.0567 


0 

0.1761 

0.2879 

0.6817 

0.2195 

0.1107 


Ca 

0.5078 

0.6638 

0.8159 

0.7184 

0.6041 

As.o. 

Ni 

0.3959 

0.9923 

1.3287 

1.3315 

0.9029 


Sn 

0.1650 

0.3300 

0.6913 

0.4757 

0.5408 


Pb 

0.2014 

0.3902 

0.6370 

0.4604 

0.4588 


DDRHF, compared with the results of RMF with PKDD [2J|. There exist significant and remarkable 
differences between DDRHF and RMF results. 

Although the attractive and repulsive parts of the nuclear force are mainly provided by a- and 
cu-mesons respectively, the contributions in DDRHF are much less than their corresponding ones in 
RMF, as shown in Fig. Q] For the isovector channels, the isovector p- and it- mesons in the DDRHF 
approach become attractive due to the stronge Fock terms. While in the standard RMF with a-, oj- 
and p-mesons, the isospin part of nuclear force is provided only by the direct part of p-meson, which 
gives the repulsive interaction for the neutrons. Furthermore it should be emphasized that one of the 
advantage of the DDRHF is the inclusion of the 7r-meson which is very important at large distance in 
DDRHF. 

An important difference between the RHF and RMF approaches is the nucleon effective mass. In 
the medium, particles or quasi-particles behave as if their mass is different from their bare mass due 
to interactions with surrounding particles, which is reflected in the level density as an example. In 
Ref. j3], the nucleon effective mass has been discussed and it is shown that there are two sources of 
modification of the bare mass: the non-locality of the mean field which gives rise to the so-called k-mass 
M7, and the energy dependence of the mean field which leads to the E-mass M E . The total effective 
mass M* is related to Mj* and M* E . One can already note that the RMF (RHF) mean field is local 
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FIG. 1: Energy density contributions from Hartree and Fock terms in different channels for neutrons in 208 Pb 
given by DDRHF with PKOl, in comparison with RMF with PKDD . 


(non-local) in coordinate space and therefore, it can be expected that their effective masses will differ. 
It should be also emphisized that, in RMF theory appears the Lorentz scalar mass Mg = M + 'Eg 
where Eg is the scalar self-energy. It should not be confused with any of the M* and one should refer 
to it as the scalar mass, or Dirac mass. 
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FIG. 2: Neutron and proton effective masses M^ R at their corresponding Fermi energy Ep calculated in 
DDRHF with PKOl as functions of /? = (N — Z)/A for different baryonic densities. 
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FIG. 3: Same as Fig.2, but for M R . 





































In the non-relativistic framework, the energy-momentum relation 


leads to the effective mass M* 



—k 2 + V(k-e) = e 
2 M y ’ ’ 


M* _ dV{k(e)-e) 
M de 


(7) 

( 8 ) 


where e = E — M is the particle energy and V(k;e ) is the momentum- and energy-dependent mean 
field. 

In a relativistic framework like RMF or RHF, the energy-momentum relation is, 


(k + kE v f + (M + S s ) 2 = (E- S 0 ) 2 , 


(9) 


where Sg, Sy, and So are respectively the scalar, spacelike- and timelike-vector components of the 
self-energy. Its Schrodinger-type form can be derived as: 


_L fc 2 + F(M _j_ 


( 10 ) 


which give the effective masses M R , 


Ml 

M 


1 de 


V ^~2M 


( 11 ) 


and = M R — e in the non-relativistic approximation by neglecting the last term at the left side 


of Eq. 


... 


One can see that M^ R is the effective mass in Refs. 



251 . |23] and M R the group mass in 


Ref. |25|], and they are the same in the non-relativistic approach. In the relativistic approach they 
can be significantly different, as shown in the following. 

The neutron and proton effective masses M^ R and M R at their corresponding Fermi energy Ep 
from the DDRHF calculations with PKOl are respectively shown as functions of (3 = (N — Z)/A 
in Fig. m and Fig. 1 for different density p At lower density, the RHF gives the trend that 
n (EF, n ) > . p (Ef,p) but this trend is reversed around 0.8po- For M R , one always have 
Mji n (EF,n) > Mp p (Ef, p ) for all densities. 

In contrast, for RMF, the relatively simple expressions tell us that one always have n < 
and M Rn (EF,n) > (Ep t p) for neutron rich system. This difference between RHF and RMF 
is related to the presence of exchange (Fock) terms which bring non-locality effects to the RHF 
self-energies. It is worthwhile to mention that in Brueckner-Hartree-Fock studies it is found that 
M*(Ep t n) > M*(Ef, p ) 3 but at larger density (pt, = 0.17fm -3 ) jroj] . 

Another significant difference between RMF and DDRHF predictions is the energy dependence 
of In RMF, Mj^p is a constant whereas it is a function of E or k in RHF. In Fig. |4j the 

energy dependence of Mj^ R for different /3 at Pb = 0.8po is shown. The neutron effective mass tends 









FIG. 4: The energy dependence of the effective mass M^ R calculated in RHF with PK01, for p = 0.8po and 
different values of neutron excess /?. 

to decrease with the energy whereas the proton mass is more constant or slightly increases in neutron 
rich matter. 

One can see that the neutron effective masses are larger than the proton ones at low energy, 
i.e., M^ Rn (E) > M^ Rp (E), and depending on the f3, a different feature appears at energy E ~ 
15 — 20 MeV ( Solid points ). It was also found in Dirac Brueckner-Hartree-Fock calculations that 
^r,„(£) > M* RtP (E) H, but at high energy (E = 50MeV). Combining with the discussion of Fig. 
|2j one can conclude that the DDRHF will predict n > M^ R at low energy or low density while 
^NR,n < -^NR,p f° r ^ le RMF. 

In summary, it has been demonstrated that one can go beyond the standard relativistic mean 
field approach to include the exchange (Fock) terms and the new couplings such as pion-nucleon 
couplings which are effective only through exchange terms. These exchange terms are the cause of 
subtle effects such as the isospin dependence of the effective masses. The same (or even better) 
quantitative description of nuclear properties comparable to RMF can be achieved with a number of 
adjusted parameters. It will open the door to the future investigation of nuclei by the relativistic 
Hartree-Fock-Bogoliubov approach and the relativistic RPA on top of RHF approximation. 
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